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Abstract 



We perform a combined analysis of the unitarity triangle and of 

the CP violating parameter e' /e using the most recent determination 

of the relevant experimental data and, whenever possible, hadronic 

matrix elements from lattice QCD. We discuss the role of the main 

,i^ . non-perturbative parameters and make a comparison with other recent 

analyses. We use lattice results for the matrix element of Qs obtained 
without reference to the strange quark mass. Since a reliable lattice 
determination of the matrix element of Qq is still missing, the theoret- 
ical predictions for e'/e suffer from large uncertainties. By evaluating 
this matrix element with the vacuum-saturation approximation, we 
typically find as central value e' /e = (4 -^ 7) x 10"''. We conclude 
that the experimental data suggest large deviation of the value of 
the matrix element of Qq from the vacuum-saturation approximation, 
possibly due to penguin contractions. 
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1 Introduction 

The latest measurements of e' je 

Re(e'/e) = (28.0±4.1) X 10-^ KTeV [1], 
Re(e'/e) = (18.5 ± 7.3) X 10-^ NA48 0], 

confirm the large value found by NA31 |^ and rise the world average to 
Re(£7e)/^^N = (21.2 ± 4.6) x 10'^ 0. Motivated by these results, we 

present a new study of the unitarity triangle and of CP violation in kaon 
decays within the Standard Model. Our results have been obtained from 
a Next-to-Leading Order (NLO) calculation of e' je combined with the con- 
straints on the Cabibbo-Kobayashi-Maskawa matrix [Vckm) derived from 
measurements of \Vcb\, \Vub\, s, AMb^ and the limits on AMb^- This work 
is an upgraded and improved version of previous studies made by the Rome 



group H]-H]. Similar analyses can be found in the recent literature []7|-|T5| 



For previous estimates of e' /e with a heavy top mass see also 116|-[20 
Several features characterize this work: 

• We analyze the constraints on Vckm together with e'/e, fully taking 
into account correlation effects. This should be compared with the 
analyses of refs. pi], E3], where only the Vckm constraints were con- 



sidered, or with the analyses of refs. [^, |T^ and |]14|, ^, where the 
input values and errors of the Vckm parameters used in the study of 
e'/e were taken elsewhere. With respect to the recent study of ref. |jl3| , 



we present the results of the analysis of the unitarity triangle together 
with the predictions for e'/e. 

In our NLO analysis, we include the full, correlated dependence of the 
coefficients of the effective Hamiltonians H^^^^, H^^^^ and H^^^^, 
computed in refs. pSf-pOf, on the relevant parameters, such as as{Mz) 



or the MS top mass, mf'^^{mf^^). 

We carefully account for the renormalization-scheme dependence of the 
matrix elements of the renormalized operators and of the compensat- 
ing effects in the corresponding Wilson coefficients. This is specially 
important, given the large differences in the values of matrix elements 



of operators defined with different prescriptions, in particular for Qq 
and Qg [} 

We address the issue of the dependence of theoretical predictions for 
e' /e on the strange quark mass m^, induced by the standard definition 
of the B parameters. In addition, we present results obtained by us- 
ing the matrix elements of the electropenguin operators, Qj and Qs, 
computed without any reference to the quark masses 



Our main results are the following. For the unitarity triangle, using experi- 
mental informations on |V^b|, \Vuh\, £, ^^3^ ^-^id AM^^, we find 

P = 0.161°:??, r7 = 0.38l°:°L 

sin(2«) = 0.1ll°:i , sin(2/3) = 0.751°;°^ , (2) 

7 = (68tl^r, JmA, = (1.34l°i^) X 10-^ 



in good agreement with the results of other recent studies |^, |22[. A detailed 
discussion of this analysis can be found in sect. 0. 

Concerning e'/e, the major uncertainty still affecting the theoretical pre- 
dictions is the lack of a quantitative determination of the matrix element of 
Qe- In particular, the status of lattice calculations for this matrix element is 
more confused now than a few years ago. The value obtained using staggered 
fermions, indicating small deviations from the vacuum-saturation approxi- 



mation (VSA) for the bare operator ^T|, has been found plagued by huge 



perturbative corrections p2] which make it unreliable. Very recent results 
using domain- wall fermions ||3^, on the contrary, correspond to a dramatic 
violation of the VSA and predict a sign of the matrix element different from 
any other non-perturbative approach. Since this result has been obtained 
with a lattice formulation for which numerical studies began very recently, 
we think that further scrutiny and confirmation from other calculations are 
needed before using it in phenomenological analyses. 

With this caveat in mind, we prefer to give firstly our result for e' /e in 
the form 

e'/e = [(-21.7lt^) GeV'' x (vrTrigf ^|ir),=o - (6.0t} j)] x 10"^ , (3) 

where the operator Q^ is renormalized at /i = 2 GeV in the 't Hooft-Veltman 
{HV) renormalization scheme of ref. [^ and the errors are evaluated by 

^ For a definition of the operators of TL^^^^^ see for example refs. po, M. 



varying all the experimental and theoretical parameters, but (7r7r|(5|^^|-ft')/=o, 
as explained in sect. |^. One may use this formula to estimate e' je with 
any non-perturbative method able to control the renormalization scale and 
scheme dependence of Qe at the NLO. 

In the absence of definite results from the lattice, we take the central value 
of l\K'n\Q^^ \K) i=Q from the VSA and allow a large variation of the matrix 
element using a relative error of 100%. Note that there is an ambiguity in 
taking this value, since the renormalization scheme and scale of the operator 
is unknown in the VSA. By assuming the VSA for the central value in the 
HV scheme at /i = 2 GeV, we find 

e'je = (3.6l|^ ± 0.5) x 10"^ (Monte Carlo), 
-11 X 10-^ < e'/e < 27 X 10"^ (Scanning), 

where the first error comes from the uncertainties on the input parameters 
and the second one accounts for the residual renormalization scheme depen- 
dence due to higher orders in the perturbative expansion. Details on the 
definition of the errors and the scanning procedure can be found in sect. ^j. 
Taking the value of the VSA in the NDR scheme, we find instead 

e'/e = {6.7tli ± 0.4) x 10"^ (Monte Carlo), 
-10 X 10^^ < e'/e < 30 X 10"^ (Scanning). 

Note that the difference between eqs. (^ and (|^) is not due to the scheme 
dependence, but to the change in the value of the matrix element of Qq. For 
a more detailed discussion of this point, see sect. |. 

We conclude that, with a central value of {717i\Qq\K)o close to the VSA 
one, even with a large error, it is difficult to reproduce the experimental 
value of e'/e, for which a conspiracy of several inputs pushing e'/e in the 
same direction is necessary. In our opinion, the important message arriving 
from the experimental results is that penguin contractions (eye diagrams), 
neglected in the VSA, give contributions to the matrix elements definitely 
larger than their factorized values. This interpretation provides a unique 
dynamical mechanism to account for both the A/ = 1/2 rule and a large 
value of e'/e within the Standard Model, whereas other arguments, as those 
based on a low value oimf'^{fi), would leave the A/ = 1/2 rule unexplained. 

This paper is organized as follows. In sect. Q, we describe the methods 
used in the phenomenological analysis. Our study of the unitarity triangle 



can be found in sect. ^ Details on the calculation of e' je are given in sect. ^. 
In both cases, we make a comparison with other recent analyses. Section ^ 
contains our conclusions. 



2 Analysis method 

The analysis of the unitarity triangle is based on a comparison of theoretical 
expressions for £, AM^^ and AM^^ with the measurements/bounds and on 
the experimental determination of |\4f,| and |V^b|. Several parameters enter 
the theoretical expressions of the above quantities. They can be classified 
in two groups: "experimental" quantities, such as the top and W masses, 
as{M^, etc., and theoretical ones, such as the Wilson coefficients, which 
are computed at the NLO in perturbation theory, and the hadronic matrix 
elements. 

We now describe the procedure used in our combined analysis of e' je and 
the unitarity triangle. 

1. For all the parameters which determine £, AM^^ and AM^^, we extract 
randomly the experimental quantities with gaussian distributions and 
the theoretical ones with fiat distributions 0. The latter include, for 
example, Bk-, the 6-quark mass, which enters as a threshold in the evo- 
lution of the Wilson coefficients, etc. We have some remarks to make 
on the choice of the error distributions. In many cases, the main sys- 
tematic error in the determination of the experimental quantities, such 
as iKibl, comes from the theoretical uncertainties. One may argue that 
systematic errors coming from the theory (and experimental-systematic 
errors as well) should correspond to fiat distributions. In practice, it 
may be difficult to disentangle statistical and systematic errors affect- 
ing some input parameters. However, as discussed below, the actual 
choice of the error distributions has a rather small influence on the final 
results. For these reasons, we have assumed for all the experimental 
quantities gaussian distributions, with a width obtained by combining 
in quadrature all errors. 

2. Among the quantities which are extracted with a gaussian distribution. 



^ In the latter case, when for a generic variable x we give the average and the error, 
X ± dx, this means that we extract x with flat probability between x — dx and x + dx. 



there are also |K 

IK- J and A = IK 



ub\ 



and iKfel; for a given set of extracted values of \Vub\, 
, we determine a = Vp^ + v'^ = |Kb|/(A|Kfe|). 



3. For a given value of a (and of all the other relevant parameters), we 
extract a value of e with gaussian distribution, and find the solutions 
of the equation 



.th. 



m 



MSf^MS 



m 



),as(M, 



Bk, a, S) 



(6) 



with respect to 6, which is the CP violation phase in the standard 
parameterization of Vckm as adopted by the PDG |^, with < S < 
TT. In eq. (|), £*''(mf^(mf^), a,(M^), . . . ,Sj^, a, 5) is the theoretical 
value computed for that given set of random parameters and e is the 
extracted value. The explicit expression of e*^ can be found for example 
in eq. (10) of ref. 0. In general one finds two independent solutions 
for 6. For any set of extractions, this fixes two independent sets of 
values for p = a cos 5 and rj = a sin 5. In the following, the set of all 
the extracted parameters and of one of the two solutions for p and rj 
will be denoted as event Q. When we consider also e'/e, event denotes 
all the random variables, including the matrix elements (and further 
parameters) which enter the calculation of this quantity. 

4. For a given event, Ui, we compute a statistical weight defined as 



Wj = exp 



1 /AM%^^-AM 



th 



rii 



2 \ dAMBa 



1 / a[AM|^^(n,)]-l 

2 [ c?a[AMf^(n,)] 



X J 



-1 



(|Kj,£*'^),(p,r/)] , (7) 



where AM^'^ (ra,) and AMj^^{ni) are the theoretical values of AMb^ 
and AMb^ computed with the event nf, a[AMB^] and da[AMBs] are 
the average value and error of the oscillation amplitude for B^-B^ 



mixing, introduced in ref. |^ f\ J MK 



'ub\ ; ^ 



th 



ip,v) 



is the Jacobian 



relating \Vub\ and e to p and rj. With this weight factor, our procedure 



^ Thus for any set of random variables we have in general two events. It may happen 
that there is no solution, in this case the event is disregarded. 

"* The values of a and da in bins of AMb^ are produced by the LEP "B Oscillation 
Working Group" . We thank A. Stocchi for providing us with these numbers. 



Constant 


Values 


Gf 


1.16639 X 10-5 GeV"^ 


ae{M,) 


7.8125 X 10-3 


sin^ 9w 


0.23154 


h 


0.1307 GeV 


Ik 


0.1598 GeV 


Mw 


80.41 GeV 


Mz 


91.1867 GeV 


Mb, 


5.2792 GeV 


Mb. 


5.3693 GeV 


Mko 


0.498 GeV 


M^ 


0.140 GeV 


AMk 


5.301 X lO^sec'i 


uo 


0.045 


RqAq 


2.7 X 10-^ GeV 


r1/2 
-D7-9 


1 


/^ 


2 GeV 



Table 1: Constants used in the numerical analysis. 



coincides with the method followed in ref. 
with flat distributions p and r/, compute 
include in the statistical weight a factor 



exp 



1 


HKr -\Vu,uv)V 


1 f e^'^'P - e{p,r])V' 


2 


^ d\Vub\ ) 


2\ de ) _ 



2^ : in that case they extract 
Kfe|(p, ??) and e{p,r]) and 



(8) 



instead of J~^. Since the error on the measurement of e is tiny, most 
of the extractions of p and t] in the interval [0, 1], correspond to very 
small statistical weights. Thus, the method of ref. ||2^ demands a 
large number of extractions in order to obtain a significative statistical 
sample. Our procedure is, in this respect, more efficient. 

5. The statistical weight Wj, suitably normalized to the sum over all the 
events, is used to compute averages and errors of the different quanti- 
ties of interest. Since the probability distributions are in general non- 
gaussian, we give the "median" and the 68% confidence level intervals. 



The median is defined in such a way that half of the weighted events 
hes below its value. The error range contains 68% of the total weighted 
events. Occasionally we will also give standard averages and errors, or 
ranges obtained by scanning the relevant parameters. A meaningful 
definition of the scanning procedure requires the "gaussian" variables 
to be extracted uniformly within some range. We choose a range ilex 
around the central values. 

We have divided the input parameters into two groups: 

A) in table |l], quantities for which the error is so small to give negligible 
effects in our analysis are listed. In this table, we also give the value of 
the renormalization scale fi at which Wilson coefficients and operator 
matrix elements are evaluated (this has not been varied) and the values 
which have been taken for B^Lg- These matrix elements have never 
been computed on the lattice and our choice is just a guess, biased by 
the VSA and justified by the fact that the precise value of these B 
parameters is not very important to the estimate of e'/e, see fig. ^j. 

B) in table |^, quantities which are extracted with gaussian (above the 
double horizontal line) or fiat (below the double line) distributions are 
listed. The value of \Vub\ is our average, and estimate of the error, of 



the CLEO |3^ and LEP |37] measurements. The lattice results comes 
from our compilation of several lattice calculations. In the table, only 
B parameters of left-left operators in the if V^-scheme are listed. For 
these operators, the physical matrix elements can be readily obtained 
by multiplying physical quantities (such as fx or M^o) times the B 
parameters, since the quark masses never enter their VSA expressions. 

Central values and errors of all the results presented in this study corre- 
spond to the "median" and the 68% confidence- level region of the appropriate 
distributions. In some cases, we also give ranges obtained by scanning all 
the input parameters. 



3 Constraints on the Unitarity Triangle 

Analyses of the unitarity triangle and NLO calculations of e'/e have been 



around for several years @]-0,[|T9|, ^, ^ ^. As discussed in the previous 



Parameter 


Value and error 


Vub 

A 

«.(M,) 

e 


0.0395 ±0.0017 

0.0037 ±0.0007 

0.2196 ±0.0023 

0.119 ±0.003 

(2.28 ±0.019) X 10-3 

(0.472 ±0.016) X lO^^sec-i 

165 ± 5 GeV 


fB^BB^/flBB, 

niB 

Bk 


4.25 ±0.15 GeV 
1.3 ±0.2 GeV 

210 ± 30 MeV 
1.14 ±0.06 
0.25 ±0.15 
0.87 ±0.13 


B parameters in the if l^-scheme at /x = 2 GeV 


Bl 
Sf 

S3 

S4 

-^5,6 
p3/2 

-07,8 


0.075 ± 0.075 

0.075 ± 0.075 

3.5 ±2.5 

3.5 ±2.5 

see eqs. (p!7|) and (|18]) 

see eq. ( 15|) 

0.63 ±0.09 



Table 2: Variable parameters: average and errors, x ± dx. The quantities 
above and below the double horizontal line have been extracted with gaussian 
and flat probability distributions, respectively. Matrix elements of left-left 
operators are expressed in terms of B parameters and measurable quantities 
(such fx, M^o, etc.) only. In these cases, it is equivalent to quote the value 
of the matrix elements or the B parameters. The latter are listed in this 
table. 



e^ 1 




Figure 1: Density plot in the p^ plane. Contours define regions containing 
5%, 68% and 95% of the generated events 



section, this analysis is based on a comparison of theoretical expressions for e, 
^Mb^ and AMb^ with the measurements/bounds and on the experimental 
determinations of \Vub\ and \Vcb\- AH relevant theoretical formulae can be 
found in refs. @-[|, 



H, HI 



Expressions, and numerical values, of the 
Wilson coefficients appearing in the effective Hamiltonians, in all popular 
renormalization schemes, can also be found in refs. P, 0, ^,p5 



__ __ The 
values of the input parameters can be found in tables |I] and |^. 

The results of our analysis of the unitarity triangle are given in eq. @. 
These values correspond to the density plot in the p-f] plane shown in fig. |I|. 

Since the fit to the unitarity triangle is overconstrained, one can extract 
one (or more) input variables, such as the renormalization group invariant 



Bk or fsa^ Bb^, together with the CKM parameters [^, ^. By removing 



the bounds on one of the above theoretical quantities, we obtain respectively 

5;, = 0.901°:?^, /b,/5^= 224+1 MeV. (9) 

These results are in very good agreement with predictions from lattice QCD, 
for example B^ = 0.89 ± 0.14 [^ (our average from a compilation of lattice 



results for Bk is given in table D and fB^yBB^ = 210±20±20 MeV || (for 



a recent review, see also ref. |0). It is worth noting that lattice predictions 
for these quantities existed long before the possibility of extracting them 
from the analysis of the unitarity triangle and have been very stable over the 
years 0. 

Comparison with other analyses 



We now briefly compare our results with the two recent analyses of refs. |2T 
and [gg. 

In spite of several differences in the procedure followed in the analysis 
of the data, our results in eq. (^ are in very good agreement with those of 
ref. [^ , which quotes Q 



P — U.iOO_oo68 ; 7 — U.O / /_o.057 ) /.^jN 

sin(2a) = OMtlH , sin(2/?) = 0.75 ± 0.090 , 7 = (671}^)'' . 



Although our method is equivalent to that of ref. [^, there are differences 
between their results, 

p = 0.202t'oml, rj = 0.340 ±0.035, 

sin(2a) = -0.26l°:l , sin(2/5) = 0.7251°;°^° , 7 = (59.51?;^)° , 



and those in eq. (0). The main reason is the range chosen for fBaV^Bd (and 
to some extent for |V"cf,|) []: in ref. p2[ the central value is similar to ours, and 



'Steve Sharpe estimated Bk = 0.84(3) (14) in 1996 |Q; a compilation by one of the 



authors of the present paper gave fBaV Bbj — 220 ± 40 MeV and fBaV Bb^ — 207 ± 30 
MeV, in 1995 and 1996 respectively g. 

^ In ref. |l|, the values of p = O.WOtomt and ry = 0.381+[5;5]^^ are given. We used 
p = p{l — A^/2) and r] — r]{l — A'^/2), with the value of A as given in table 0. 

^ We thank A. Stocchi for pointing it out to us. 
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Figure 2: Density plot showing the correlation between p and /b^v-Bb, 



to the value used by ref. ^1 
larger values of /^^ 



B 



B, 



but the error is strongly asymmetric favouring 
This choice is justified by arguing that the existing 



unquenched results for fs^ are larger than the corresponding quenched ones. 
This is not a good argument, however, because the effect of quenching on 

the B parameter, i.e. on the full matrix element parametrized by fs^yBB^, 
is still unknown and could compensate the increase of /s^. Thus, for this 
quantity, we have chosen the symmetric range in table |^. Given the strong 

correlation between p and /s^y^s^, as shown in fig. ^ this choice accounts 
for the bulk of the differences in the results. We checked that, with a similar 

range for fs^y Bb^, our results are much closer to those of ref. ^2 . 



Besides other minor differences in the range of the input parameters, the 
various analyses also differ in the treatment of the errors. The authors of 



ref. 1^ attempt to distinguish, for each input parameter, statistical and sys- 
tematic errors and assign to them different distributions. In ref. [21|, all the 
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relevant input variables are included in the expression of the x^-function to 
be minimized, implicitly assuming gaussian distributions for all these quan- 
tities. In our case, we have defined two classes of statistical- and systematic- 
dominated parameters (see table 0), and use a single distribution, either 



gaussian or fiat, for each variable. Finally, both refs. [^Jl and ||2^ consider 
the Wilson coefficients relevant to K and B meson mixing as independent 
input parameters, assuming for them either gaussian or flat distributions. 
This is not appropriate since the QCD corrections are known functions of cts, 
mt, etc. and not independent quantities. In our analysis, for each event, we 
have computed the NLO QCD corrections, thus taking correctly into account 
the correlations. 



Apart from the choice of the range of fB^yBB^, these further differences 
have small effect on the final results for the unitarity triangle. In addition, 



the fitted values of Bk and fBaV Bb^ in eq. (§) are well consistent with the 



1-0.34 



results of ref. [g2|, where they found Bk = 0.871^;^^ and Jb^^Bb^^ = 223 ±13 
MeV (Scenario I), and with those of ref. [^, I3k = O.SOlgjg and fB^yBB^ = 
222^11 MeV. 

4 Calculation of e' /e 



In this section, we summarize the main steps necessary to the calculation 
of e'/e and give some details on the procedure followed to obtain the result 
quoted in eq. (PP-(^. We also make a critical comparison with previous 



calculations of the Rome group [H-[H, with the recent estimates of ref. |13 



and with results obtained in other approaches such as the xQM B and 



the 1/N expansion [|T^, pO 



Some general remarks are necessary before entering a more detailed dis- 
cussion. Given the large numerical cancellations which may occur in the 
theoretical expression of e'/e, a solid prediction should avoid the "Harlequin 
procedure" . This procedure consists in patching together Bq from the xQMy 
Bg from the 1/N expansion, m^'^{fi) from the lattice, etc., or any other 
combination/average of different methods. All these quantities are indeed 
strongly correlated (for example Bq and Bg in the 1/A^ expansion or B pa- 
rameters and quark masses in the lattice approach) and should be consis- 
tently computed within each given theoretical framework. Unfortunately, 
as will clearly appear from the discussion below, no one of the actual non- 
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perturbative methods is in the position to avoid completely the Harlequin 
procedure, not even for the most important input parameters only. The 
second important issue is the consistency of the renormalization procedure 
adopted in the perturbative calculation of the Wilson coefficients and in the 
non-perturbative computation of the operator matrix elements. This prob- 
lem is particularly serious for the xQM and the 1/A^ expansion, and will 
be discussed when comparing the lattice approach to these methods. We 
will address, in particular, the problem of the quadratic divergences appear- 
ing in the 1/A^ expansion. This is an important issue, since the authors of 
refs. |T^, 1^ find that these divergences provide the enhancement necessary 
to explain the large values of Re^o and of e' je. 
Schematically, e' can be cast in the form 



, exp(i7r/4) uj 
V2 ReAo 



cj-^ImAa - (1 - fi/ij)ImAo (12) 



where lo = ReA2/Rey4o and ReAg, given in table |T|, are taken from exper- 
iments, and QjB is a correcting factor, estimated in refs. |^-|^, due to 
isospin-breaking effects. Using the operator product expansion, the K -^ tctt 
amplitudes lraA2 and ImAo are computed from the matrix elements of the 
effective Hamiltonian, expressed in terms of Wilson coefficients and renor- 
malized operators 

{n7r\n^''^='\K') = -^Y.C^{^^){^^\Q^{^^)\K') (13) 

where the sum is over a complete set of operators, which depend on the 
renormalization scale //. In general, there are ten four-fermion operators and 
two dimension-five operators representing the chromo- and electro-magnetic 
dipole interactions. In the Standard Model, the contribution of these dimension- 
five operators is usually neglected (possible SUSY effects can enhance the 
contribution of the chromomagnetic operator [Q) and, with the scale /i = 2 
GeV > rr^c at which calculations are performed ^ only 9 out of the 10 four- 
fermion operators are independent p6|, |2^. Wilson coefficients and matrix 
elements of the operators Qi{fi), appearing in the effective Hamiltonian, sep- 
arately depend on the choice of the renormalization scale and scheme. This 
dependence cancels in physical quantities, such as Im742 and ImAo, up to 



^ /i = 2 -;- 3 GeV > rric is the typical scale at which matrix elements are computed in 
lattice QCD. 

13 



1 o 
O.l 



-0.04 — 




Q,° Q/ Qj Q. Qs Qs Qy'^ Q„'^^ Q,"^ Qy""^ Q/'''' Q. 



Figure 3: Individual contribution of the different operators to e' je, using the 
value of {Qq)q in eq. (17). The corresponding value ofe'/e is also shown. 



higher-order corrections in the perturbative expansion of the Wilson coef- 
ficients. For this crucial cancellation to take place, the non-perturbative 
method used to compute hadronic matrix elements must allow a definition of 
the renormalized operators consistent with the scheme used in the calculation 
of the Wilson coefficients. 

So far, lattice QCD is the only non-perturbative approach in which both 
the scale and scheme dependence can be consistently accounted for, using 
either lattice perturbation theory or non-perturbative renormalization tech- 
|50| . This is the main reason why we have followed this approach 



49 



niques 

over the years. 

In fig. H, we display the individual contributions of each operator to the 
prediction ofe'/e using matrix elements computed on the lattice, as explained 
in the next section. Qualitatively, very similar results are obtained also in 



14 



other approaches WU, O, M. There is a general consensus that the largest 



contributions to e'/e are those coming from Qq and Qs, with opposite sign, 
and sizeable contributions may come from Q4, Q^ and Qq in the presence of 
large cancellations between Qq and Q%, i.e. when the prediction for e'/e is 
of (9(10~^) 0. For this reason the following discussion, and the comparison 
with other calculations, will be focused on the determination, and errors, of 
the matrix elements of the two most important operators. In the following, 
we adopt the following notation: 

(Q,){o,2} = (^7r|Qf^(/i = 2 GeV)|ir°),={o,2} (14) 

where the superscript denotes the HV renormalization scheme as defined in 
ref. [H Q. 

Status of the calculation with matrix elements from lattice QCD 

The evaluation of physical K ^ tttt matrix elements on the lattice relies 
on the use of Chiral Perturbation Theory (xPT): so far only {TT\Qi{fi)\K) 
and {t^{p = 0)7r(g = Q)\Qi{ii)\K)i=2 (with the two pions at rest) have been 
computed for a variety of operators. The physical matrix elements are then 
obtained by using xPT at the lowest order. This is a consequence of the dif- 
ficulties in extracting physical multiparticle amplitudes in Euclidean space- 



time ||5T[. Proposals to overcome this problem have been presented, at the 
price of introducing some model dependence in the lattice results [^. The 
use of xPT implies that large systematic errors may occur in the presence of 
large corrections from higher-order terms in the chiral expansion and/or from 
final-state interactions (FSI). This problem is common to all approaches: if 
large higher-order terms in the chiral expansion are indeed present and im- 
portant, any method claiming to have these systematic errors under control 
must be able to reproduce the FSI phases ^o and 82 of the physical ampli- 
tudes. The approaches of ref. |^, [10[ and |]n|, however, give FSI smaller 



than their physical values. Regarding this issue, we note that the idea of 
improving the predictions of the hadronic amplitudes using the experimental 
values of the FSI phases, with formulae such as {Qi)i — *■ Re{Qi)j/ cos 6'j^^, 
is illusory. If the FSI phases are not theoretically under control, one cannot 



^We recall the reader that Qs and Qg give large contributions only to the 1=2 amplitude. 
^° Incidentally, we note that the HV scheme of ref. p5|] is not the same as the HV 



scheme of ref. 
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tell whether the main uncertainty conies from the real or the imaginary part 
of the computed amplitude, or from the absolute value needed to compute 
Aq and A2. 

1 = 2 matrix element of Qs- There exists a large set of quenched calcula- 
tions of {Q8)2 performed with different formulations of the lattice fermion ac- 
tions (Staggered, Wilson, tree-level improved, tadpole improved) and renor- 
malization techniques (perturbative, boosted perturbative, non-perturba- 
tive), at several values of the inverse lattice spacing a~^ = 2 -h 3 GeV Il3( 

>(3/2) 



o^,ra-55 



All these calculations, usually expressed in terms of Bg ' , give 
consistent results within 20% of uncertainty. Among the results, we have 
taken the central value from the recent calculation of ref. |^, where the 
matrix elements {Q8)2 and ((57)2 have been computed directly without any 
reference to the quark masses, and inflated the errors to account from the 
uncertainty due to the quenched approximation (unquenched results are ex- 
pected very soon) and the lack of extrapolation to zero lattice spacing. The 
values we use are 



(g7)2 = 0.18 ± 0.06 GeV^ 
(Qg)2 = 0.62 ± 0.12 GeV^ 



(15) 



The operator matrix elements computed without reference to quark masses 
are given in physical units. The reader who likes to work with B parameters 
and quark masses may use the formulae 



(Q8(/i))2 



V2U 
V2U 



1 



Ml, 



Mlo - Ml 



3 Vmf^(/i)+mf^(^), 



Ml, 



Ml, - Ml 



,mf'5(/i) +mf^{^), 



Br\^), 



The values of the matrix elements in eq. ( |15|) correspond to B^ = 0.89 ± 
0.30 and ^f ^^^ = 0.93±0.18 for a "conventional" mass fixed to mf^+mf^ = 
130 MeV at /i = 2 GeV. Anybody may rescale the value of the B parameters 
for these operators according to her/his preferred value for m^^^{ii). 

Matrix element of Qq. For {Qq)q from the lattice, the situation appears 
worse today than a few years ago when the calculations of refs. [§[-0 were 
performed: 
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i) until 1997, the only existing lattice result, obtained with staggered 
fermions (SF) without NLO lattice perturbative corrections, was Bq = 
1.0 ± 0.2 [pT|. This is the value used in our previous analyses M, R|; 



ii) with SF even more accurate results have been quoted recently, namely 
Bq = 0.67 ± 0.04 ± 0.05 (quenched) and B^ = 0.76 ± 0.03 ± 0.05 (with 
Uf = 2) S; 

iii) 0{as) corrections, necessary to match lattice operators to continuum 
ones at the NLO, are so huge for Qq with SF (in the neighborhood of 
— 100% 1^) as to make all the above results unreliable. Note, however, 
that the corrections tend to diminish the value of (Q6)o; 

iv) the latest lattice results for this matrix element, computed with domain- 
wall fermions from {7i\Qq\K) |]5^, are absolutely surprising: {Qg)o has 
the sign opposite to what expected in the VSA, and to what is found 
with the xQM and the 1/A^ expansion. Moreover, the absolute value is 
so large as to give e'/e ~ —120 x 10~^. Were this confirmed, even the 
conservative statement by Andrzej Buras [^, namely ... that certain 
features present in the Standard Model are confirmed by the experimen- 
tal results. Indeed the sign and the order of magnitude of e'/e predicted 
by the Standard Model turn out to agree with the data... would re- 
sult too optimistic. In order to reproduce the experimental number, 
e'/e ~ 20 X lO"'^, not only new physics is required, but a large can- 
cellation should also occur between the Standard Model and the new 
physics contributions. Since this result has been obtained with domain- 
wall fermions, a lattice formulation for which numerical studies began 
very recently, and no details on the renormalization and subtraction 
procedure have been given, we consider it premature to use the value 
of the matrix element of ref. [Q in phenomenological analyses. Hope- 



fully, new lattice calculations will clarify this fundamental issue. 

v) estimates in the framework of the 1/A^ expansion (where, however, one 
should always to take into account the correlation between the values 
of Be and sf ^^^ [||]) and by the xQM ^ give Bq = 0.7^ 1.3 at scale 
/i = 0.6 -f- 1.0 GeV. One may argue fl^, ^ that the scale dependence 
of Bq above 1 GeV is rather weak and take the range Bq = 0.7 -r- 1.3 
as valid also at /i = 2 GeV, which is the scale at which we work. 
Note, however, that the dependence of the matrix elements on the 
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Density plots showing the correlation of Im Xf and {Qq)o with 



renorinalization scheme is rather strong and that, in these approaches, 
the scheme in which matrix elements are computed is unknown. 

Taking into account i)-v), we conclude that there is no computation of (Qq) 
that can be reliably used in phenomenological analyses. For this reason, we 
present our result as in eq. @. In addition, biased by the VSA, we give results 
assuming B^^ = 1.0 ± 1.0 and B^^ = 1.35 ± 1.35 (which corresponds to 
Bq^^ = 1.0 ± 1.0), taking an error of 100%. We compute the corresponding 
matrix element with the same "conventional" quark mass used for {Qs)2, for 
which an explicit lattice calculation without reference to the quark masses 
exists. In physical units, this choice corresponds to 



(Qe), 



-0.4±0.4GeV^ 



(17) 



and 

(gg)o = -0.6±0.6GeV^ (18) 

and (^5)0 = l/3(Q6)o) in the two cases. These expressions are obtained 
using the formulae 

4 / /\/f2 \^ 

(Q5(^))0 = -- K^=— ^^=— (/x-/.)55(/i), 



3 VmfS(/i)+mf^(/i), 



2 



<*<"»» ^ -' L^s^jAtf^j ^''^ - ^-'^''^^ ■ <^'' 



The values of e' /e obtained using the matrix elements in eqs. (|T7|) and (18) 
are given in eqs. (^) and (^) respectively. The large difference between the 
two results is due to the strong correlation between e' /e and (Qe); as shown in 
fig. ^. In the same figure we also show the correlation between e' /e and ImXt, 
which is the parameter governing the strength of CP violation in i^-meson 
decays. 

In fig. ^ we give two event distributions of e'/e, obtained using Wil- 
son coefficients and hadronic matrix elements consistently computed in HV 
and NDR. Indeed, the calculation of the Wilson coefficients achieved in 
refs. [^, 0, 1^, [^-[0 allows a consistent determination of the matrix ele- 



ments of the renormalized operators at the NLO. Notice that the two distri- 
butions in fig. 1^ are quite similar, while differences between matrix elements 
in different schemes can be rather large: the parameter Bq decreases, for in- 
stance, by ~ 30% going from the HV to the NDR scheme. This decrease is, 
however, largely compensated by a readjustment of the corresponding Wil- 
son coefficients, as must happen in NLO calculations. The uncertainties due 
to higher order perturbative corrections, given by the second error in eqs. (^ 
and (^, have been evaluated by modifying consistently Wilson coefficients 
and matrix elements in the HV and NDR schemes. In the two cases, using 



for example (Q6)o from eq. ([T7|), we obtain 



e'/e = (3.1+IDxlO-^ HV , 

e'/e = (4.0+^i) X 10"^ NDR, (20) 

from which the result in eq. (^) has been derived. 

The results in eqs. (^ and (|^) are in very good agreement with previous 
estimates of the Rome @l-0 and Munich group 0, |^. This agreement it 
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Figure 5: Probability distribution of e' /e in the HV and NDR schemes. In 
changing schemes matrix elements and Wilson coefficients have been consis- 
tently redefined at the NLO. 



is not surprising since the two groups are using very similar inputs for the 
matrix elements and the experimental parameters have only slightly changed 
in the last few years. The crucial problem, namely a quantitative determina- 
tion of {Q6)o, remains unfortunately still unsolved. At present, we can only 
conclude that, with the central value of {Qg)o taken from the VSA, even 
with a large error, it is difficult to reproduce the experimental value of e'/e. 
On the other hand, by scanning various input parameters {as{Mz), ImXt, 
etc. and, in the conventional approach, Bq and Sg ) and in particular by 
choosing them close to their extreme values, it is possible to obtain e'/e up to 
30 X lO""*. This gives the impression of a better agreement (lesser disagree- 
ment) between the theoretical predictions and the data. For completeness, 
we also give the interval of values of e'/e obtained by scanning, within one a 
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the different parameters. We obtain 



- 11 X 10"^ < e' je < 27x 10-\ 

-10 X 10"^ < e'je < 30 X 10"^ , (21) 



using {QojQ from eqs. ( |17|) and ([18|) respectively. Equation (l2l|) allows a 



direct comparison with several calculations appeared in the literature 

0.2 X 10"^ < e' je < 22.0 x 10"^ Munich99-HV [13], 
1.1 X 10-^ < e' je < 28.8 x 10"^ Munich99-NDR, 
7.0 X 10"^ < e' je < 31 x 10"^ Trieste98 



(22) 



1.5 X 10-^ < e' je < 31.6 x 10^^ Dortmund99 |4 



Note that our scanning results include a region of negative e'/e on account 
of our choice of the error on (Q6)o being larger than in other cases. 

In spite of the fact that the experimental world average is compatible 
with the above "scanned" ranges, we stress that, in order to get a large 
value of e'/e, a conspiracy of several inputs pushing e'/e in the same di- 
rection is necessary. For central values of the parameters, the predictions 
are, in general, much lower than the experimental results. For example, the 
Rome-Munich and 1/A^ estimates are typically in the range 3-8 x 10~^ and 
8-10 X 10~^, respectively. For this reason, barring the possibility of new 
physics effects [HSl, we believe that an important message is arriving from 



the experimental results: 

penguin contractions (or eye diagrams, not to be confused with penguin op- 
erators ffU^), usually neglected within factorization, give contributions to the 
matrix elements definitely larger than their factorized values. 

This implies that the "effective" B parameters of the relevant operators, 
specifically those relative to the matrix elements of Qi and Q2 for ReAo and 
of Qe for ^'/^ ^^^ much larger than one. This interpretation would provide a 
unique dynamical mechanism to explain both the A/ = 1/2 rule and a large 
value of e'/e [^ . Large contributions from penguin contractions are actually 
found by calculations performed in the framework of the Chiral Quark Model 
(xQM) [0, [g or the l/N expansion |T3|, |2g, §21-|g. It is very important 



that these indications find quantitative confirmation in other approaches, for 
example in lattice QCD calculations. Note that naive explanations of the 
large value of e'/e, such as a very low value of mf^^{fi), would leave the 
A J = 1/2 rule unexplained. 
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We have quantified the amount of enhancement required for the matrix 
element of Qg in order to explain the experimental value of e' je. A fit of 
(<56)o to the world average Re{e' le)wA-, using for the other parameters the 
standard values given in tables |l] and|| (in particular by varying {Q%)2 in the 
interval given by eq. (|15|) , since for this operators penguin contractions are 
absent), gives (Q6)o = — l-2lo;ii ± 0.15 GeV^, about 2^3 times larger than 
the central values used in our analysis. In units which are more familiar to the 
reader, the value of (Q6)o required to fit the data corresponds to -Be = 2 -i- 3 
for mf^ + m^ = 105 -i- 130 MeV (in the HV scheme). 

Before ending this discussion, we wish to illustrate the correlation existing 
between the B parameters and the quark masses in lattice calculations. 

On the lattice, quark masses are often extracted from the matrix elements 
of the (renormalized) axial current (A^) and pseudoscalar density {P{fi)) (for 
simplicity we assume degenerate quark masses) 

^ 1 {a\d,A,\P) 
'"^^^-2(«|P(/i)|/5)' ^^^^ 

where a and /3 are physical states (typically a is the vacuum state and (3 
the one-pseudoscalar meson state) and m(/i) and PifJ^) are renormalized in 
the same scheme. On the other hand, the B parameters of Qq and Qs are 
obtained (schematically) from the ratio of the following matrix elements, 
evaluated using suitable ratios of correlation functions |^: 

where P-jr and Pk are the pseudoscalar densities with the flavour content of 
the pion or kaon, respectively. Eqs. ( p3[ ) and (p^ demonstrate the strong 
correlation existing between B parameters and quark masses: large values of 
the matrix elements of -P(/u) correspond, at the same time, to small values 
of m(/i) and i?6,8(A^)- Physical amplitudes, instead, behave as 

(Q6,8) = const. X ^^ , (25) 

where "const." is a constant which may be expressed in terms of measur- 
able quantities (specifically Mk and fx) only. From eqs. (^3]) and (|2^) , we 

^^ See for example ref. [S^l- We omit the superscript (3/2) in B^ for simplicity. 
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recognize that the dependence on {P{fi)) cancels in the ratio BQ^s/m{fiY, 
appearing in the physical matrix elements. 

Previous lattice studies preferred to work with B parameters because 
these are dimensionless quantities, not affected by the uncertainty due to the 
calibration of the lattice spacing. This method can still be used, provided 
that quark masses and the B parameters from the same simulation are pre- 
sented together (alternatively one can give directly the ratio BQ^^/m{iif'). 
In ref. |^0|, two possible definitions of dimensionless "5 parameters", which 



can be directly related to physical matrix elements without using the quark 
masses have been proposed. In this analysis we have used the values of 
{Q7,s)i=2 computed with one of these new definitions. 

Comparison with the Munich group 

The original approach of the Munich group was to extract the values of 
the relevant matrix elements from experimental measurements |^, §. This 
method guarantees the consistency of the operator matrix elements with the 
corresponding Wilson coefficients. In this approach, a convenient choice of 



the renormalization scale is /x = nic- In their analysis the authors of ref. ||T^ 
used the value fi = nic = 1.3 GeV. 

Unfortunately, with the Munich method it is impossible to get the two 
most important contributions, namely those corresponding to {Qg)o and 



{Q8)2- In this respect, we completely agree with ref. |T3| that one cannot 



extract {Qe)o from the experimental value of ReA^, unless further assump- 
tions are made [0. For this reason, "guided by the results presented above 
and biased to some extent by the results from the large-A^ approach and 
lattice calculations", the authors of ref. [^ have taken Bq = 1.0 ± 0.3 and 



(3/2) 

Bg = 0.8 ± 0.2, at /i = 1.3 GeV. These values, if assumed to hold in the 
HV regularization, are close to ours, given the smooth behaviour of the B 
parameters between /i = 1.3 and 2.0 GeV. The main differences in the eval- 
uation of {Qe)o and {Q8)2 between ref. [^ and our calculation come from 
the value of mf^i^n) and from the scheme dependence that we now discuss. 
In a complete NLO calculation, the scheme dependence of the matrix 
elements is compensated by that of the Wilson coefficients up to NNLO 
terms, so that physical quantities are independent of the renormalization 
procedure. Lattice QCD allows a complete control of the definition of the 
renormalized operators at the NLO: for example, {Q%)2 has been computed 
with specific NLO definitions. In ref. [|13|, however, they kept fixed the 
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values of the B parameters when changing the renormahzation scheme of the 
Wilson coefficients. Although it is true that {Qe)o is, at present, unknown, 
this procedure introduces an unphysical scale and scheme dependence which 
should be avoided. The use of the same B parameters in two different schemes 
leads to overestimate the error due to the scheme dependence. We consider 
more appropriate to increase the error on the matrix element {B parameter) 
in a given scheme and attribute the final uncertainty to our ignorance on the 
matrix element, rather than to the choice of the renormahzation scheme. 

For comparison, we also present the results obtained by taking the main 
parameters close to those of ref. |T^, namely B^^^ = 1.0±0.3, -Bg — 



0.8 ±0.2, mf^{^i) = (110 ±20) MeV and the condition Bq > Bf''^\ We find 

e'/e = (7.2+1^) x 10"^ , (26) 

a central value well consistent with the result of ref. [|l^ in the NDR scheme, 
given the remaining differences in the renormahzation scale and in the other 
matrix elements. However, we were not able to reproduce the long posi- 
tive tail in the e'/e distribution of ref. [|1^], which produces an error more 



asymmetric than that in eq. (^6]) . 



1/A^ expansion and xQM 

The 1/A^ expansion and the xQM are effective low energy theories which 
describe the hadronic world. To be specific, in the framework of the 1/N ex- 
pansion the starting point is given by the chiral Lagrangian for pseudoscalar 
mesons expanded in powers of masses and momenta. At the leading order in 
1/A^ local four-fermion operators can be written as products of currents and 
densities, which are expressed in terms of the fields and coupling of the effec- 
tive theory. In higher orders, in order to compute the relevant loop diagrams, 
a (hard) cutoff, Ac, must be introduced. This cutoff must be lower than ~ 1 
GeV, since the effective theory only includes pseudoscalar bosons and can- 
not account for vector mesons or heavier excitations. The cutoff is usually 
identified with the scale at which the short- distance Wilson coefficients must 
be evaluated. 

Divergences appearing in factorizable contributions can be reabsorbed 
in the renormalized coupling of the effective theory. Non-factorizable cor- 
rections constitute the part which should be matched to the short distance 
coefficients. By using the intermediate colour-singlet boson method, the au- 



thors of refs. |T^, |6^ claim to be able to perform a consistent matching, 

24 



including the finite terms, of tlie matrix elements of the operators in the 
effective theory to the corresponding Wilson coefficients. It is precisely this 
point which, in our opinion, has never been demonstrated in a convincing 
way. 

If the matching is "consistent" , then it should be possible to show that 
in principle the cutoff dependence of the matrix elements computed in the 
1/A^ expansion cancels that of the Wilson coefficients, at least at the order 
in 1/A^ at which they are working. Moreover, if really finite terms are under 
control, it should be possible to tell whether the coefficients should be taken 
in HV, NDR or any other renormalization scheme. 

The fact that in higher orders even quadratic divergences appear, with 
the result that the logarithmic divergences depend now on the regulariza- 
tion, makes the matching even more problematic. Theoretically, we cannot 
imagine any mechanism to cancel the cutoff dependence of the physical ampli- 
tude in the presence of quadratic divergences, which should, in our opinion, 
disappear in any reasonable version of the effective theory. Note that, in 
refs. 1^, |T^ , the calculations are performed using dimensional regularization 
in which quadratic divergences do not appear. We suggest to repeat the cal- 
culation of the relevant matrix elements with the xQM using a hard cutoff to 
show the stability of the results with respect to the change of regularization 
scheme and verify the possible presence of quadratic divergences. This would 
also provide an easier comparison with the 1/A^ expansion. 

It is also necessary to show (and to our knowledge it has never been 
done) that the numerical results for the matrix elements are stable with 
respect to the choice of the ultraviolet cutoff. This would also clarify the 
issue of the routing of the momenta in divergent integrals. For example, 
the matrix elements in the meson theory could be computed in some lattice 
regularization. 

5 Conclusions 

In this paper, we have presented a combined analysis of the unitarity trian- 
gle and e'/e, using (whenever possible) matrix elements from lattice QCD, 
which is theoretically well suited for NLO calculations. We stress that, given 
the correlations among different non-perturbative parameters, calculations 
of e'/e should use matrix elements consistently computed within a given the- 
oretical approach, at least for the main contributions. At present, however, 
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there is no reliable calculation of {Qq)o on the lattice. For this reason, our 
main result is the one given in eq. (j^) , in which the matrix element is left as a 
free parameter. In addition, we give the results in eqs. (^) and (|^), by taking 
the central value of {Qq)o from the VSA in different renormalization schemes 
with a relative error of 100%. Our results show that, even with such a large 
error, it is difficult, although not impossible, to reproduce the experimental 
value of e' /e. To this end, a conspiracy of several input parameters, pushing 
in the same direction, is necessary. We rather think that the important mes- 
sage arriving from the experimental results is that penguin contractions (eye 
diagrams) give contributions which make the matrix element of Qq definitely 
larger than what expected on the basis of the VSA. This interpretation pro- 
vides a unique dynamical mechanism to account for both the A/ = 1/2 rule 
and a large value of e' /e within the Standard Model, whereas other argu- 
ments, as those based on a low value of m^^^ , would leave the A/ = 1/2 
rule unexplained. Concerning the strange quark mass, we stress that m^^ is 
irrelevant for the calculation of the operator matrix elements on the lattice. 
In the long run, lattice QCD is the only non-perturbative method able to 
produce quantitative results at the NLO accuracy. In the present situation, 
however, other approaches, such as the xQM or the 1/A^ expansion, which 
cannot control the proper definition of the renormalized operators, may prove 
useful to understand the underlying dynamics. We hope that the issue of the 
computation of {Qi,2) and {Qe)o on the lattice will be clarified soon, providing 
reliable theoretical estimates for both the A/ = 1/2 rule and e'/e. 
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